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INTRODUCTION 


Kith  roughly  ten  thousand  events  of  m.  >  4.0 
occurring  each  year,  the  possioility  exists  that  short- 
period  signals  from  an  explosion  or  an  earthquake  will 
be  masked  in  the  coda  of  another  earthquake.  When  the 
signals  are  masked  at  only  a  few  stations,  detection 
is  often  possible.  If,  however,  the  short-period  sig¬ 
nals  from  an  event  are  masked  at  all  stations,  it  seems 
likely  (due  to  the  high  background  level  on  long-period 
records,  and  to  the  persistence  of  motion  from  large 
events)  that  the  long-period  signals  will  also  be 
masked.  If  this  occurs,  detection  is  not  possible. 
Obviously,  then,  some  attempt  should  be  made  to  search 
the  secondary  portions  of  seismograms  for  phases  which 
may  indicate  the  presence  of  a  new  event. 

While  simple  in  principle,  the  search  of  the  coda 
of  an  event  for  secondary  arrivals  is  fraught  with 
difficulties.  Though  Cohen,  et  al .  (1972)  have  demon¬ 
strated  that  coda  characteristics  are  controlled  pri¬ 
marily  by  the  arrival  times  of  predictable  secondary 
phases  (c.g.,  ppf  PcP,  etc.),  numerous  signal-like 
excursions  within  a  given  coda  cannot  be  explained  in 
terms  of  least-time  paths  through  the  earlh.  These 
unexplained  phases  may  be  thought  of  as  "alarms"  in 
the  sense  that  their  existence  will  cause  a  more 
thorough  search  of  the  pertinent  seismograms  to  be 
made  for  other  indications  of  a  masked  event.  Most 
often  the  phase  will  go  unexplained.  In  a  larger 
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context ,  however ,  we  would  almost  certainly  want  evi¬ 
dence  of  a  second  c^cnt  on  the  seismograms  at  four  or 
more  stations,  so  that  the  event  epicenter  can  he 
determined.  Ihc  ability  to  locate  the  second  event  is 
important  if  coda  suppression  techniques  such  as  the 
beam  and  mixed-signal  processor  (Dean  ct  al.,  1968) 
arc  to  be  used  for  enhancement  of  the  masked  signal, 
ihus,  we  arc  concerned  with  location  solutions  obtained 
using  random,  unexplained  phases.  Specifically,  we 
arc  concerned  with  the  number  of  significant  location 
determinations  obtained  for  nonexistent  events;  that 
is,  the  false  alarm  probabilities  for  mixed  events. 

In  determini  ,g  false  alarm  probabilities,  we  pro¬ 
ceeded  as  follows.  About  1500  P  and  PKP-codas  were 
searched  for  phases  which  exceed  a  given  amplitude 
threshold.  These  phases  were  first  purged  of  known 
secondary  arrivals  (c.g,  PP,  PCP,  etc.).  The  ratio 
of  the  number  of  unexplained  arrivals  to  the  number 
of  coda  examined,  then,  determined  the  probability  of 
an  unexplained  phase  occurring  at  a  single  station. 

The  binomial  theorem  was  subsequently  used  to  deter¬ 
mine  the  probability  P4  of  having  unexplained  arri¬ 
vals  in  a  given  event’s  coda  at  four  stations  out  of 
a  world-wide  network  of  N  stations.  That  these  unex¬ 
plained  arrivals  yield  a  significant  location  solution 
will  occur  with  probability  Pu<  The  product  of  P4 
and  PL4  then  yields  the  probability  that  four  stations 
of  the  world-wide  network  will  experience  a  false  alarm. 
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METHOD 


Assume  that  the  probability  of  an  unexplained 
phase  occurring  in  a  coda  of  length  T  is  PQ.  Then, 
for  a  network  of  N  stations,  the  probability  P  that 
four  or  more  stations  will  have  unexplained  phases 
for  a  given  event  is: 


p5  .  P  <(,P  )N, 


5(1-P  )N"5 
o  v  o 


(1) 


■  P4 


♦ 


♦ 


•  •  • 


If  PLN<N-4>  is  thc  probability  that  unexplained 
phases  observed  for  a  given  event  at  N  stations  will 
yield  a  significant  location  solution,  thc  probability 
PFA  a  ^a^se  alarm  is 

PFA  *  P4-PL4  *  P5*PL5  *••• 
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DATA 


ihc  f^rst  quantity  to  determine  is  the  probability 
PQ  ot  an  unexplained  phase  occurring  at  a  single  sta¬ 
tion.  To  do  this,  we  use  the  P  and  PKP  coda  determina¬ 
tions  compiled  by  Cohen,  et  al.  (1972).  Specifically, 
for  "anivals"  which  exceed  a  given  threshold  ampli¬ 
tude,  we  will  examine  the  number  of  unexplained  phases 
present  in  the  coda. 

The  method  used  to  determine  the  coda  decay 
characteristics  is  shown  in  Figure  1.  Amplitude  mea¬ 
surements,  sealed  relative  to  the  largest  excursion 
in  the  P  or  PKP  coda,  were  made  in  a  specified  set  of 
successive  time  windows.  Measurements  were  made  until 
the  coda  decayed  into  the  pre-existing  ambient  noise 
level,  or  until  a  period  of  fifteen  minutes  had  elapsed 
from  the  time  of  cbda  onset.  In  a  few  eases,  measure¬ 
ments  terminated  with  the  arrival  of  a  second  event. 

The  principal  coda  maxima  were  next  plotted  on  log- 
linear  paper,  and  the  coda  envelope  obtained  by  con¬ 
necting  successive  determinations.  For  example,  the 
coda  measurements  of  Figure  1  yield  the  coda  envelope 
shown  in  Figure  2.  Using  the  records  from  seventeen 
World  Wide  Standard  Seismograph  Stations  whose  loca¬ 
tions  are  shown  in  Figure  3,  Cohen,  et  al.  (1972) 
determined  the  coda  characteristics  for  events  from 
the  fifteen  seismic  regions  shown  in  Figure  4j  examples 
of  typical  coda  decay  characteristics  arc  shown  in 
Figure  5. 


The  following  example  illustrates  the  method 
used  to  flag  suspect  phases  in  the  coda.  As  shown  in 
Figure  2,  we  assume  the  log-coda  decays  linearly,  and 
draw  a  line  to  connect  the  observations  plotted  at 
40  and  60  seconds  elapsed  time.  The  relative  ground 
motion  at  50  seconds  is  then  predicted,  yielding  a 
value  of  about  30%.  The  observed  value,  however,  is 
50%;  as  such,  the  signal-to-coda  ratio  r  (observed-to- 
predicted)  is  1.67  (4.5  db) .  Had  we  set  a  detection 
threshold  ratio  r  of  1.5  (3.5  db) ,  the  arrival  plotted 
at  SO  seconds  would  be  flagged  for  subsequent  analysis. 
If  this  arrival  is  not  found  tc  oe  a  predictable 
secondary  phase  (i.e.,  pP,  PP,  PcP,  SKP,  etc.),  it  is 
considered  a  false  alarm.  Obviously,  recourse ’to  the 
travel- time  tables  for  body  wave  and  surface  wave 
Phases,  and  to  the  actual  seismograms,  is  required  to 
purge  the  flagged  arrivals  of  predictable  phases. 

Note  that  no  relationship  exists  between  the  detection 
threshold  r  and  relative  event  size.  To  see  this  we 
need  only  consider  that  the  earthquake  and  the  event 
masked  may  be  at  different  distances  from  the  receiver, 
or  that  the  event  masked  may  appear  anywhere  in  the 
coda  of  the  earthquake. 

Examples  of  false  alarm  analyses  are  shown  in 
I ables  I  and  II.  By  computing  all  observed-to- 
predicted  coda  ratios,  cumulative  figures  for  unex¬ 
plained  phases  may  easily  be  obtained  for  any  given 
threshold  ratio.  In  this  work  we  analyze  only  those 
phases  for  which  the  detection  ratio  is  >  1.5.  Using 
1,471  P  and  PKP  coda  determinations,  coda  search  and 
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analysis  techniques  yield  the  unexplained  phases 
listed  in  Tables  III  through  XVII.  In  general,  these 
phases  can  be  ascribed  to  the  following  phenomena: 

A.  Pulsing  (Figure  6a) 

Pulsing  is  taken  to  mean  any  large  (in  comparison 
to  the  background  level)  phase  or  excursion  in  the 
coda  which  is  not  a  predictable  arrival. 

B.  Normal  Coda  decay  characteristics  (Figure  6b) 

Coda  amplitudes  rise  and  fall  in  the  natural 
course  of  decay.  These  undulations  arc  often  large 
enough  to  produce  an  unexplained  phase  as  we  have  de¬ 
fined  this  term.  An  example  of  this  phenomena  is  the 
set  of  coda  determinations  (70t,  701,  301),  which 
yields  a  detection  ratio  of  1.60  (4.1  db) .  Natural 
decay  characteristics  in  a  coda  arc  generally  identi¬ 
fied  as  such  by  an  experienced  analyst,  and  would 
probably  be  dismissed  from  further  consideration. 

To  be  uniform  in  our  treatment  of  the  data,  we 
include  all  unexplained  phases  -  regardless  of  their 
character  -  in  our  compilation. 

A  summary  of  the  phase  analyses  is  shown  in  Table 
XVIII.  For  a  detection  threshold  ratio  of  1.5  (3.5  db) , 
the  data  yield  175  unexplained  phases  in  1,471  events. 
Thus,  the  probability  PQ  of  such  a  phase  occurring  at 
a  single  station  is  0.118  or  roughly  one  unexplained 
phase  for  every  ten  events.  If  we  consider  only  those 
detections  for  which  the  threshold  is  2.0  (6.0  db) , 


the  single-station  probability  drops  to  .00816,  or 

less  than  one  unexplained  phase  for  every  100  events. 

In  both  cases  the  average  coda  length  T  is  343  seconds 
(-6  minutes) . 


Taolc  XX  summarizes  the  results  for  the  random 
location  analyses.  For  any  given  trial,  a  random  set 
of  four  stations  is  chosen  from  the  available  world¬ 
wide  network  (complete  n« twork  shown  in  Table  XIX). 
For  each  set  of  four  stations,  random  times  between  0 
and  T  (T  -  6  or  10  minutes)  are  added  to  the  arrival 
times  associated  with  an  event  at  the  center  of  the 
region  specified.  The  altered  arrival  times  are  then 
used  to  determine  a  new  surface- focus  epicenter.  If 
the  absolute  value  of  each  station's  residual  travd- 
time  error  |ei|  is  less  than  a  specified  value  (here, 
or  6  seconds),  the  location  is  considered  acceptable, 
"significant”.  As  an  example,  the  analysis  performed 
for  the  Tadzhik-Hindu  Kush  region  using  10  minute  coda 
yields  50  acceptable  solutions  in  4505  trials  (|e  |  < 

6  seconds);  thus,  Pu  -  0.011.  The  significant  1  " 

location  solutions  are  shown  in  Figure  7;  by  chance, 
one  of  the  location  solutions  is  for  a  non-existent 
event  in  Novaya  Zemlya.  The  randomly  derive^  location 
mobt  remote  from  the  reference  epicenter  (34’ ;j,  73-c) 

1  located  5582  km  to  the  northwest,  at  coordinates 
( S 5 °N ,  6°E) . 


3 

or 
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RESULTS 


Results  of  the  false  alarm  analyses  are  summarized 
in  Table  XXI.  In  computing  falsc-alarm  probabilities, 
we  assume  that  a  thirteen-station  network  is  used  for 
monitoring  purposes  (sec  Appendix). 

lor  a  detection  threshold  of  r  «  1.5  (3.5  db) ,  a 

6  minute  coda,  and  a  value  l^l  6  seconds,  a  thirteen- 

station  world-wide  network  would  experience  1.5  false 
alarms  per  1,000  events.  With  roughly  10,000  events 

of  magnitude  m^  4.0  occurring  each  year,  the  network 
would  experience  15  false  alarms  per  year.  However, 
many  of  the  events  with  magnitudes  on  the  order  of  4.0 
will  have  coda  lengths  much  less  than  6  minutes,  re¬ 
sulting  in  fewer  false  alarms.  If  jej  <  3  seconds, 
the  false  alarm  rati  dre"-  by  a  factor  of  2,  to  about 

7  false  alarms  per  year. 

False  alarm  rates  for  10  minute  coda  (  r  ■  1.5) 
arc  similai  to  the  above.  While  the  probability  PL4 
of  four,  random  unexplained  phases  yielding  a  signifi¬ 
cant  location  solution  decreases  with  increasing  coda 
length,  the  assumption  or  a  uniform  distribution  in 
the  occurrence  of  unexplained  phases  results  in  an  in¬ 
crease  in  P4.  Thus,  the  data  of  Tab  1»-  XXI  suggest 
that  for  r  *  1.5  and  an  increase  in  roda  lengths  to 
10  minutes,  the  product  of  the  probabilities  P4  and 
PL4  incrcascs  slightly,  but  by  a  factor  less  than  2.0. 
For  10  minute  coda  and  JeJ  <  6  seconds,  then,  we 
expect  about  17  false  alarms  per  year.  For  |e^|<  3 


seconds,  the  false  alarm  rate  drops  to  9  false  alarms 
per  year. 

At  a  detection  threshold  of  2.0  (6  db) ,  false 
alarms  arc  virtually  eliminated. 

While  the  probabilities  P^4  determined  for  the 
Tadzhik-Hindu  Kush  region  were  used  for  the  computa¬ 
tions  of  false  alarm  piobabilitics  (Table  XXI),  use 
of  the  PL4  determined  for  the  Philippine  Islands  or 
the  Kamchatka  Peninsula  would  not  significantly  alter 
the  results. 

The  false  alarm  probability  for  four  or  more 
stations  is  essentially  the  same  as  that  for  exactly 
four  stations.  To  sec  this,  consider  a  case  for  coda 
lengths  of  6  minutes  and  r  -  1.5;  from  Table  XXI,  we 
note  that  PQ  -  0.118.  Thus,  for  a  thirteen  station 
worldwide  network,  we  find  from  equation  (i)  that 

Ps  “  0,24  P4*  Now»  requiring  that  five  stations  with 
random  arrival  times  yield  a  significant  location  solu 
tion  is  a  more  severe  restriction  than  requiring  four 
stations  to  produce  such  a  location.  Thus,  as  sec^tn 
Tables  XXa  and  b  for  the  Kamchatka  Peninsula  dati* 

PL5  =  *04  P4  (lcil  1  6  seconds).  Taken  together,  the 
probability  that  a  four  or  five  station  network  will 
experience  a  false  alarm  is  approximately  II  greater 
than  the  probability  of  a  four  station  network  experi¬ 
encing  a  false  alarm.  Contributions  for  sub-networks 
of  six  or  more  stations  will  be  even  less.  Thus,  our 
estimates  for  the  probability  of  a  false  alarm  at  four 


stations  are  for  all  practical  purposes  the  same  as 
the  falsc-alarm  probability  for  four  or  more  stations. 
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DISCUSSION 


Basic  to  the  problem  of  determining  false  alarm 
probabilities  for  mixed  events  is  the  question  of  what 
detection  threshold  is  to  be  used  for  phase  detection, 
he  have  seen  that  for  a  coda  length  of  6  minutes,  a 
threshold  of  r  *  1.5  (3.5  db)  yields  about  15  false 
alarms  per  year  for  IcJ  £  6  seconds  and  7  false  alarms 
per  year  for  lej  <  3  seconds.  At  r  -  2.0  (6.0  db) , 
false  alarms  are  virtually  eliminated.  It  appears, 
then,  that  a  network  can  be  operated  with  r  ■  1.5,  or 
perhaps  slightly  lower,  without  unreasonable  false 
alarm  occurrence. 

It  should  also  be  noted  that  the  number  of  unex¬ 
plained  phases  occurring  in  the  P  coda  at  the  various 
stations  of  a  world-wide  network  is  quite  sensitive  to 
the  choice  of  a  threshold  level.  As  we  have  seen  from 
Table  XVI 1 1 ,  decreasing  the  threshold  from  r  -  2.0  to 
r  *  1.5  increases  the  number  of  unexplained  phases  by 
more  than  an  order  of  magnitude  (12  versus  175  unex¬ 
plained  phases  for  1,471  coda  observations).  These 
observations  suggest  that  the  threshold  used  for  phase 
detection  in  the  coda  of  an  event  be  in  the  neighbor¬ 
hood  of  r  ■=  1.5,  but  probably  no  lower. 

Given  that  time  picks  made  in  the  secondary  por¬ 
tions  of  seismograms  are  subject  to  considerable  eiror, 
even  a  consistent  set  of  arrival  times  for  a  masked 
event  may  yield  an  estimated  location  error  of  several 
liundred  kilometers.  For  reasons  related  to  on-site 


inspection,  therefore,  let  us  inquire  into  the  loca¬ 
tion  uncertainty  implied  by  residual  travel-time  errors 
ei  for  which  |c.|  <6  and  3  seconds. 

In  determining  the  area  associated  with  specified 
arrival- time  errors,  we  randomly  perturb  the  known 
travel  times  for  a  given  event  to  four  stations  selec¬ 
ted  at  random  from  a  world-wide  network,  and  examine 
the  scatter  in  the  location  solutions.  Details  of 
this  analysis  arc  given  in  the  Appendix.  From  the 
Appendix,  we  find  the  travel-time  errors  6.  for  which 
"6  1  6i  1  6  imply  that  an  event  can  be  located  any- 
uherc  in  an  area  3.2  x  105  kn.2  in  size  (95%  confidence 
level).  Further,  the  location  analysis  shown  in  the 
Appendix  suggests  that  errors  in  the  travel-time  picks 
less  than  or  equal  to  | 6 i J  imply  residual  travel-time 
eirors  ^  such  that  |e.|  :  |6.|/2.  Thus,  the  3.2  x  105km2 
area  of  unccrtainity  should  be  associated  with 
residual  travel-time  errors  -3  <  e.  <3  seconds.  If 
I c i I  -  6  seconds,  | 6 i |  <  12  seconds,  and  the  area  of 
unccrtainity  is  about  4  times  larger,  or  1.3  x  106  km2. 
Similarly,  if  |c.|  <  1.5  seconds,  |6.|  <  3  seconds, 
and  the  area  of  unccrtainity  is  roughly  8.0  x  104  km2. 
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CONCLUSIONS 


1.  Using  1,471  p  an(j  pkp  coda  determinations  (Cohen, 
et  cl • ,  19/2)  coda  analysis  /iclds  the  following: 

Detection  Phases 

threshold  (r)  flagged 

1.5  (3.5  db)  990 

2  0  (6.0  db)  446 

In  general,  unexplained  phases  resulted  from  pulsing 
(unpredictable  excursions  in  the  coda)  and  normal  coda 
undulations.  The  average  coda  length  for  the  events 
studied  was  about  6  minutes  (543  seconds). 

2.  The  probability  Pq  of  an  unexplained  phase  occur¬ 
ring  in  a  6-minute  coda  at  a  single  station  is  0.12 
for  r  ■  1.5  (3.5  db).  Thus,  the  probability  P4  that 
four  stations  out  of  thirteen  will  have  unexplained 
phases  in  the  coda  of  a  given  earthquake  is  0.045.  As 
the  probability  that  unexplained  phases  at  four  stations 
v  1.  yield  a  significant  location  solution  if  0.032 

for  residual  travel-time  errors  -3  £  c.  ^3  seconds, 

the  probability  of  a  false  alarm  occurring  is  0.0015 

per  event.  This  is  equivalent  to  about  15  false  alarms 
each  year. 

3.  The  false  alarm  probability  for  four  or  more  sta¬ 
tions  is  essentially  the  same  as  the  false  alarm  prob¬ 
ability  for  four  stations. 

4.  For  a  detection  threshold  of  2.0  (6.0  db) ,  false 
alarms  are  virtually  eliminated. 


Unexplained 

Phases 

175 

12 
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>.  The  threshold  used  for  phase  detection  should  be 
in  the  neighborhood  of  r  *  1  5  (3.5  db) ,  but  probably 
no  lower. 

(» .  Tor  purposes  of  on-site  inspection,  residual 
travel-time  errors  "3  £  3  seconds  imply  that  an 

event  can  be  located  r  /where  in  an  area  3.2  x  105  km2 
in  size  (95%  confidence  level). 

.  If  residual  travel-time  errors  are  such  that 
-1.5  £  £  1.5  seconds,  the  false  alarm  rate  drops  by 

a  factor  of  2,  to  7  false  alarms  per  year  for  r  *  1.5. 
The  corresponding  area  of  uncertainity  for  location 
purposes  is  reduced  by  a  factor  of  4,  to  8.0  x  104  km2 
(951  confidence  level). 
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figure  3.  Map  showing  location  of  worldwide  network. 
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gurc  4.  Regions  used  in  coda  ana 


n«f 

»vr. 

<10  Ofv 

0*< 

0.0 

0 

♦  *. 

r  •  7* 

r. 

10.0 

s 

IT. 

0.  1  7 

i". 

70.0 

0 

H. 

0.11 

i". 

10.0 

0 

47. 

0.11 

i". 

♦  0.0 

01. 

o.|0 

10. 

00.0 

*4. 

o.l  1 

K. 

1.0 

*0. 

0.1* 

1". 

1.0 

♦  1. 

".1* 

|0. 

/ 1  * 

70. 

".IT 

*. 

"•* 

71. 

0.71 

T. 

*  0  J 

7T. 

0.70 

T. 

70. 

0.7* 

». 

®  •  # 

7*. 

0.  io 

0. 

»•  * 

11. 

0.1* 

0. 

^1  * 

7T. 

0.71 

0. 

*•  * 

IT. 

0.74 

0. 

10.0 

70. 

0.1* 

♦. 

t  1.0 

10. 

o.io 

♦. 

17.0 

II. 

0.01 

7. 

11.0 

10. 

0.7* 

7. 

1  • 

•  • 


11  *•  ••  II  Ht  III  HI  l|| 

TIMI  If  III  IMIVII  (MIIHMtl 

Figure  5a.  P-coda,  Alaska,  CMC. 
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Figure  5.  Coda  characteristics. 
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Significant  location  solutions, 
nute  coda,  -6  <  c.  <  o  seconds) 


Phase  Analysis,  P-Coda,  CMC,  Alaska, 
11  April  1966  (18:26:11.82) ,  A  -  20° 


TIML  INTI  RVAL 
s:  seconds 
m:  minutes 

RATIO 

(Observed-To- Predicted) 

COMMENTS 

5-  10s 

0.79 

10- 20s 

1.34 

20- 30s 

0.92 

30-40s 

1.25 

40- 50s 

0.58 

50-60s 

1.71 

Unexplained  phase 

1  -  2m 

0.82 

2-  3m 

1.01 

3-4n 

0.75 

4  -  5m 

1.63 

SS ,SSS , PeP, Surface  Waves 
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tabu;  ii 

PBasc  Analysis,  l>KI'-Coda,  WES,  Philippine 
Taiwan,  21  November  1965  (10:31:49.72),  A 


Islands  - 
-  139° 


T i ml  interval  ratio 

s:  seconds  (Observcd-to-Predicted) 
m:  minutes 


COMMENTS 


5-  10s 

0.98 

10- 20s 

0.80 

20- 30s 

1.15 

30-40s 

0.57 

40- 50s 

1.88 

50-60s 

0. 71 

1  -  2m 

1.19 

2-  3m 

0.59 

3-4m 

1.83 

4-5m 

0.92 

5-6m 

0.78 

6-  7m 

1.43 

pPKP 


SKP ,  PP 
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TABLE  III 

Unexplained  Phases  -  South  America 
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a. 
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TABLH  V 

Unexplained  phases  -  California  and  Western  United  States 
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TABLL  VI 

Unexplained  Phases  -  Alaska 


TABLL  VI 1 

Unexplained  Phases  -  Aleutian  Islands 
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APPENDIX 

Computation  of  the  Average 
Location  Area  Associated 
With  Specified  Arrival  Time  Errors 
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i'c  seek  a  determination  of  the  area  associated 
uitn  specified  errors  in  station  arrival  times. 

Fo‘  „  given  region  outlined  in  Figure  4  (main 
text)  the  distances  from  the  four  corners  of  the  region 
to  each  station  listed  in  Table  XIX  (main  text)  were 
computed.  If  the  distance  between  any  one  corner  and 
a  station  exceeds  100°.  this  station  is  eliminated 
irom  the  network.  The  P  wave  travel  times  from  the 
center  of  the  region  to  the  remaining  stations  in  the 
network  are  then  computed.  Using  a  random-number 
generator,  four  stations  arc  selected  from  the  network, 
and  their  travel  times  perturbed  by  random  numbers  of 
up  to  ±6.  seconds.  Program  SHIFT  (Chiburis,  1968)  is 
then  used  to  obtain  a  surface-constrained  location. 

In  all,  20  location  solutions  using  random  sets  of 
four  stations  and  random  travel-time  perturbations  are 
computed  for  each  region. 


Of  these,  one  is  deleted,  yielding  an  estimate 
(at  the  95*  confidence  level)  of  the  location  area 


Aj..  Averaging  over  the 
mate  J  for  the  location 


^k  yields  a  world-wide  csti- 
arca.  A  typical  plot  of  the 


locations  obtained  for  random  combinations  of  four 


stations  and  random  sets  of  time  perturbations  is 


given  in  Figure  A1  (Solomon  Islands  -  New  Hebrides). 
Here,  we  have  assumed  that  arrivals  can  be  picked  to 


an  accuracy  of  better  than  6  seconds,  that  is, 

^  -  ^i  1  ^  °ads .  We  delete  the  solution  at  (6.5°S, 
152.8  L),  and  enclose  the  remaining  locations  with  a 


rectangle  centered  at  the  true  epicentral  location 
1  ^ * S »  156°E)*  This  rectangle  defines  the  area  A.  in 
which  95t  of  all  lorjtions  are  assumed  to  fall.  The 
area  in  Figure  A1  is  7‘  high  by  3‘  wide.  Applying 
short-distance  conversion  factors  (Richter,  Appendix 
UI,  1958),  the  area  is  found  to  be  775  km  high  by 


325  km  wide;  the  area  A^,  therefore,  is  2.5  x  10! 


km' 


A  summary  of  the  solutions  for  the  Ak  is  given  in 
Table  AI.  As  noted,  an  insufficient  number  of  stations 
was  available  for  analysis  of  the  South  American  region, 
further,  the  determinations  for  the  Philippine  Is.- 
Taiwan  region,  and  Turkey-Greece  are  questionable,  as 
less  than  19  solutions  were  available  for  analysis. 

Using  the  results  tor  the  12  remaining  regions  yields 
an  average  area  *  of  3.2  x  105  km2,  with  an  average 
of  13  stations  receiving  P-wave  arrivals  from  each 
region  (-6  <  6.  <  6  seconds).  For  -3  <  c.  <  3  seconds, 
the  area  dimensions  given  in  Table  AI  should  be 
roughly  halved.  This  may  be  seen  from  Table  All  which 
shows  the  results  of  two  sets  of  location  analyses 
performed  using  random  travel-time  perturbations  of 
up  to  ±3  and  ±6  seconds.  For  a  given  region,  the 
same  20  sub-networks  of  four  stations  were  used,  but 
the  travel-time  perturbations  in  each  trial  for  which 
-3  <  Ci  <  3  seconds  were  exactly  half  of  those  used 
in  the  corresponding  trial  for  which  -6  <  6.  <  6  seconds. 
In  each  case,  the  area  dimensions  for  -3  <  6.  <3  seconds 
are  about  half  of  the  dimensions  obtained  for 
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-0  £  6.  6  seconds.  Thus,  tor  picking  errors  of  up 

to  t3  seconds,  the  average  location  area  X  is  on  the 

c  2 

order  of  0.8  x  10  km  . 

As  a  final  note,  when  determining  the  location 
of  an  event,  the  easiest  way  to  estimate  the  accuracies 
of  the  travel-time  picks  would  seem  to  be  by  examina¬ 
tion  of  the  residual  travel-time  errors  for  each 
station.  To  determine  the  experience  relationship 
between  residual  travel-time  errors  and  picking 
errors  6^,  let  us  compare  these  parameters  for  various 
location  solutions  computed  using  *-pccificd,  random 
travel-time  perturbations. 

Figure  A2  shows  a  comparison  of  residual  travel- 
time  errors  and  the  specified  (input)  picking 
errors  6^.  Kc  examine  19  location  solutions  for  each 
of  two  regions;  in  both  eases,  ”6  <_  6.  <  6  seconds. 

Of  the  152  data  pairs  shown,  142,  or  93%  of  the  total, 
lie  within  the  area  bounded  by  -3  <  <  3  seconds. 

A  similar  study  for  which  -3  <.  6^  3  seconds  (Figure 

A3)  shows  about  901  of  the  determinations  to  fall  in 
the  range  -1.5  1.5  seconds.  Thus,  the  data 

suggest  that  travel-time  errors  less  than  or  equal  to 
|6^|  imply  residual  travel-time  errors  c.  such  that 
Ujl  =  l«il/2. 
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*  Omit  from  computation  of  averages 

**  obtained  from  data  far  which  -3  6  *3  seconds  and  adjusted  to  -6  <  £  <  6  seconds. 


tabu;  aii 


Corresponding  Location  Analysis 
for  Pickinj  l.rrors  4.  of  up  to  13  and  16 


Seconds 


Pickirg  Errors 

Up  to  ±3  seconds 
Up  to  ±6  seconds 

Up  to  i3  seconds 
Up  to  ±6  seconds 


Area  Dimensions* 
Confidence  Level) 
Lat i tude  Longitude 


1.0°  1.4° 

1- 4°  2, 8® 

2-  1°  2.9° 

4.7°  5.9° 


Region 

Kamchatka  Peninsula 
Kamchatka  Peninsula 

Kurile  Islands 
Kurile  Islands 


20  solutions  computed;  one  deleted 


r 
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lOOiliVI  Ninos 


EAST  LONGITUDE  |i«|) 


figure  Al.  Location  solutions, 
Hebrides  (12°S,  156°Ii). 


Solomon  Islands 


New 


-  A  8  ■ 


figure  A3.  Comparison  of  residual  travel-time  errors  and 
arrival-time  picking  errors  (-3  <  6^  <  3  seconds). 


